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'I'he  GtatuG  of  the  current  work  cu  ruby  Iriuer  induced  fan  ionisation 
is  reported.  iixperimental  reacure'r.entc  of  the  focused  spot  dia;'eter, 
breakdown  cnission,  ionization  tirr.es  and  nonattcriuation  ionization  ar»- 
jiresented.  Tiie  theoretical  study  results  obtained  to  date  in  tiie  areas 
of  nonlinear  ionization  initiation  and  the  extension  of  the  rderc'wave 
breakdown  theory  are  discussed. 
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I.  ■[NTRODnCTTni! 


k.  PURPOSE 

This  is  the  second  technical  documentar;'’  report  for  contract  mmVer 
AF30(  ('Ol  )333.?  under  the  sponsorship  of  the  Rome  Air  Development  Tenter.  he 
period  of  work  covered  is  from  11  August  19^i+  to  11  December  191*.  The  status 
of  the  current  work,  results  obtained  to  date,  and  future  •■/ork  to  be  performed 
are  discussed. 

■'his  contract  is  for  the  investigation  and  study  of  the  breakdovm  (ioniza¬ 
tion)  mechanisms  of  gases  at  optical/infrared  freo/iencies.  The  investigation 
includes  both  theoretical  and  experimental  research.,  using  lasers,  into  the 
power  densities  required  to  ionize  air  betvreen  abc'.t  1000  pounds  per  3q''are  inch 
an  !  1C  ^  nillir.eters  of  mercury,  and  the  mechanism, s  responsillc  for  c  :ch  ioniza¬ 
tion.  Approaches  employed  acco'unt  for  mechanisms  of  breakdown  of  several  p  ire 
'•ases  and  gas  mixt’ires  chosen  y  the  researchers,  as  vrell  as  air.  Ionization 
is  determined  as  a  f'unction  of  tine,  gas  ionization  potentials,  pressire,  fre¬ 
quency  of  incident  beam,  presence  of  outside  agents  (container  -.Alls,  surfaces, 
ionizing  radiation),  and  other  pertinent  parameters. 

B.  PROGRAM 

\s  shovm  in  the  gas  ionization  program  flow  diagram,  figure  I-l, 
parallf’l  programs  vrre  conducted,  one  for  experL~ents  and  the  other  for  theoreti¬ 
cal  studies.  A  continuous  interchange  exists  betvm'en  the  experimental  and 
study  areas,  permitting  theory  to  help  guide  the  experimental  approach  and  the 
experimental  results  to  help  develop  and  confirm  the  theoretical  descriptions. 
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Figure  I-l.  Gas  Ionization  Program 


’’asiwUy,  two  types  of  experiments  are  to  be  performed:  (1)  the  measure¬ 
ment  of  the  net  ionization  rates  and  (2)  the  statistical  time  lags.  The 
critical  nature  of  the  experiments  required  first  an  experimental  setup 
calibration.  For  program  '•xpediency,  the  net  ionization  rates  experiments 
are  divided  into  tvro  phases.  The  first  phase  determines  those  rates  for  a 
few  representative  gases  at  selected  gas  pressures  and  allows  the  basic 
mechanism  to  be  isolated  early  in  the  program,  thus  aiding  the  theoretical 
analysis.  The  second  phase  extends  these  rate  measurements  for  the  full  range 
of  press'ires  and  gases. 

"he  first  sot  of  net  ionization  rate  experiments  js  to  determine  the  rate 
of  ionization  buildup  once  the  process  is  initiated.  If  it  is  exponen¬ 
tial,  the  nechanipm  can  be  characterised  as  an  avalanche-type  process 
similar  to  the  well-known  microwave  breakdown.  However,  if  the  buildup  is 
linear,  the  ionization  rate  m\jst  depend  on  a  different  type  of  mechanism.  At 
tb'"  high  power  densities  required  for  optical  breakdown  very  high  field 
strong*  .3  exist.  These  fields  nay  be  sufficient  to  induce  direct  ionization 
of  the  gas  molecules  by  photons,  rather  than  depending  upon  the  collision  mechanism 
that  accounts  for  microwave  breakdown. 

The  statistical  time  lag  experimants  are  to  determine  hov  the  ionization 
process  is  initiated.  The  avalanche  ionization  buildup,  which  depends  upon 
initial  ions  being  present  within  the  volume  subject  to  the  electric  field, 
furnishes  the  electrons  that  start  the  ionization  process.  The  random  nature 
of  the  ion  distribution  in  the  gases  should  cause  a  statistical  time  lag  in  the 
Icnization  initiation  for  those  mechanisms  that  depend  upon  their  presence.  If 
the  intense  optical  field  can  directly  ionize  the  gas  molecules,  a  much  less 
statistical  effect  in  the  ionization  initiation  will  be  observed. 


The  theoretical  studies  are  divided  into  three  areas.  Tvo  of  the  study 
areas  bet^an  in  the  early  part  of  the  program  (i.e.,  the  Investigation  of  the 
non-linear  ionization  initiation  mechanisms  and  the  development  and  extension 
of  the  ’  reakdovm  theory).  The  third  study  area  correlates  the  expexdjnental 
res’ilts  with  breakdown  theory.  If  disagreement  is  evident,  the  theory  will  be 
modified  to  agree  with  the  experiments  to  obtain  a  quantitative  description  of 
the  breakdown  mechanisms. 

r..  3"'ATUS 

\t  the  end  of  this  reporting  period,  measurements  of  the  focused  spot 
diameter  for  several  lenses  were  completed.  Data  on  breakdown  emission  from 
Argon  during  and  after  laser  Irradiation  wore  obtained.  Ionization  times 
versus  electric  field  for  different  focal  spot  diameters  were  also  recorded  for 
Argon.  Preakdown  measurements  were  obtained  for  halltou  Charge  collection 
measurements  were  made  at  very  low  air  pressures  to  try  to  determine  the 
ionization  build-up  process.  licperiinents  have  begun  using  a  variable  pressure 
test  coll. 

The  theoretical  studies  considered  both  the  non-linear  ionization  initiation 
mechanisms  and  the  extension  of  the  microwave  breakdown  theory.  Attempts  to 
find  useful  expressions  for  the  rates  of  non-linear  ionization  have  continued 
but  witliout  success.  Investigation  of  the  rate  of  photon  absorption  in  free-free 
transitions  has  continued. 
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II.  3M-IWIY  Airo  C0NCi;j3I0N3 


A.  F:0CE35  DESCRIPTION 

Tlie  theoretical  and  ojcperimental  data  available  at  present  are  consistent 
with  the  following  sequence  of  events.  It  appears  that  the  growth  of  ionization 
begins  with  the  production  of  a  small  number  of  electron-ion  pairs  from  LTpuri- 
ties  or  from  the  main  gas  by  a  non-linear  or  nulti-photon  process.  The  electron 
and  ion  density  then  begins  to  grow  by  an  electron  multiplication,  or  avalanche 
process,  which  probably  involves  raultiphoton  photoionization  of  excited  atoms. 

■,'hen  a  very  small  fraction  of  the  gas  is  ionized  (e.g.,  less  than  one  part  in 
10  at  atmospheric  pressure),  space  charge  effects  become  large  enough  to  reduce 
t(.e  electron  loss  and  increase  the  rate  of  electron  multiplication,  Firvally, 
w^ien  the  gas  is  nearly  con^iletely  ionized,  the  absorption  of  laser  energy  in 
electron-neutral  collisions  is  augmented  by  absorption  in  electron-ion 
collisions  and  in  some  cases  by  absorption  by  excited  atoms.  This  final  stage 
produces  the  observed  absorption  of  the  laser  beam  and  most  of  the  light  output. 

,  .  pm  DENTAL  WORK 

The  process  of  gas  ionization  shows  a  threshold  dependence  upon  the  in¬ 
tensity  distribution  in  the  laser  beam.  Since  significant  changes  in  the  in¬ 
tensity  distribution  can  cause  significant  changes  in  spot  size,  care  must  be 
taken  to  monitor  these  variations.  Secondly,  spot  size  for  short  focal  length 
lenses  cannot  be  detemined  from  a  knowledge  of  beam  divergence,  thus  making 
accurate  measurements  of  spot  size  an  iJz^)o^tant  parameter  to  define  experimentally. 
Reducing  the  spot  size  by  use  of  a  5.08  cm  focal  length  lens  lowered  the  time 
of  breakdown  for  a  given  field  strength  as  con5)ared  to  that  for  a  3.3  cm  lens. 
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results  of  the  vrork  done  In  helium  Indicate  that  the  spectral  emission 
of  the  plasma  differs  for  the  various  gases.  The  easily  ionizable  mercury 
vapor  had  no  detectable  effect  on  the  breakdown  time  in  helium. 

'T!’',nR"'”IC\L  ’.'ORK 

Evidence  is  presented  that  the  absorption  of  relatively  high  energy 
photons  in  free-free  transitions  may  in  some  cases  lead  to  sra-ill  (A/ 30"^)  Vmt 
possibly  significant  deviations  in  the  rate  of  growth  of  an  electron  avalanche 
from  that  predicted  by  an  appropriate  extension  of  the  microwave  theory,  i.e., 
the  continuous  form  of  the  Boltzmann  equation.  This  evidence  is  obtained  by 
solving  the  difference  form  of  the  Boltzmann  equation  in  both  the  steady  state 
and  time  dependent  cases  for  various  simple  models  of  the  electron-atom  or 
electron-ion  collision  processes,  liowever,  since  direction  and  magnitude  of 
the  deviation  depend  upon  the  model  used  and  since  several  iir^xjrtant  effects 
are  not  included,  further  conqMirisons  of  the  two  theories  are  desirable. 

Additional  calculations  of  the  excitatior.^  ionization,  and  diffusion 
coefficients  are  presented  and  used  to  compute  avalanche  growth  times  and 
approximate  curves  of  the  laser  intensity  required  for  breakdown  as  a  function 
of  gas  density.  In  general  the  agreement  between  theory  and  experiment  is  within 
a  factor  of  two  in  laser  intensity  or  40^  in  laser  electric  field  intensity. 

There  are,  however,  systematic  deviations  between  theory  and  experiment  which 
should  be  explored  further. 
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III.  r,K3  rONIZ/.TION  kxpi=:rih?iits* 


VAHIA-’inNS  LI  THPL''.H0LD  KKAStmLMTJlTrj 

The  measurenentr?  of  the  time  of  breakdovm  In  argon  as  determined  from 
the  time  of  attenuation  of  the  transmitted  laser  light  versus  total  energy 
were  shown  in  Figures  III-9,  IIT-10,  and  III-U  of  the  first  interim  report. 

In  those  measurements  a  large  difference  In  the  threshold  of  breakdown  (lowest 
energy  producing  attenuation)  existed  between  Figures  111-9  and  III-IO  11. 
Differences  in  the  data  due  to  two  separate  gas  cells  were  ruled  out  Ijy 
interchanging  the  cell  used  to  picoduce  the  data  shown  in  Figure  III-9  and 
observing  data  almilar  to  that  shown  in  III-IO  and  III-U.  Replacement  of 
the  focusing  lens  was  then  thought  to  have  caused  the  discrepancy.  Further 
analysis  of  the  intensity  distribution  showed  the  real  cause  of  the  threshold 
change.  An  analysis  of  the  intensity  distribution  had  shown  that  the  effective 
beam  divergence  varied  between  2.1  and  2,6  milliradians .  At  that  time  it  was 
thought  that  this  beam  divergence  was  characteristic  of  the  system  at  all 
times.  The  analysis  of  the  be«un  divergence  during  later  ejqperlmentation  revealed 
that  the  angle  of  divergence  was  then  about  0.8  milliradians.  These  changes 
in  beam  divergence  were  observed  during  subsequent  expezdjaents  and  measured  to 
vary  by  as  much  as  several  milliradians.  The  changes  seemed  to  accompany  the 
replacement  of  the  xenon  flashlanq)  used  to  excite  the  ruby.  The  variation  in 
the  angle  of  divergence  during  the  life  of  any  one  lanqp  seemed  to  be  about  25%. 

As  yet  there  does  not  seem  to  be  a  noticeable  correlation  between  these  smaller 
variations  in  divergence  angle  and  the  snail  variation  in  time  of  breakdown 
versus  total  energy. 

*Thl8  work  was  perfonaed  by  R.  W.  Waynant 
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l-EASURIMKMT  OF  FOCUSED  SPOT  SIZE 


The  measurement  of  the  focused  spot  size  was  undejrway  at  the  time  of  the 
first  Interim  report.  The  measuring  technique  was  outlined  in  that  document. 

The  data  for  several  lenses  has  since  been  analyzed  and  will  be  reported  here. 

Graphs  showing  the  diameter  of  the  spot  as  a  function  of  distance  from  the 
center  of  the  lens  are  shown  in  Figures  III-l,  III-2,  and  III-3.  In  the  case 
of  the  lenses  with  12.7  cm  and  5.08  cm  focal  lengths  a  glass  plate  was  inserted 
between  the  lens  and  the  photographic  film  in  order  to  simulate  the  conditions 
encountered  with  the  gas  cell.  This  probably  accounts  for  the  fact  that  the 
smallest  measured  diameter  was  foimd  at  distances  greater  than  the  focal  length 
stated  by  the  manufacturer.  Since  the  3.3  cm  focal  length  lens  will  go  into 
the  stainless  steel  test  cell,  the  spot  size  for  two  of  these  lenses  were  measured. 

In  all  cases,  the  msasnred  dlhaiter  of  the  spot  ilia  was  larger  than  the  diameter  cal¬ 
culated  from  angular  dlrergenee  data.  This  indicates  that  spherical  abberatlon 
in  the  lens  plays  an  important  part  in  spot  size  determinations. 

The  measurements  of  spot  size  for  the  12.7  cm  lens  show  considerable 
varlat{on  of  the  spot  diameter  versus  distance  from  the  lens.  These  variations 
as  shown  in  the  plot  have  not  been  corrected  for  changes  in  the  intensity  dis¬ 
tributions,  but  such  correction  seems  only  to  exagerate  the  variations.  This 
would  perhaps  suggest  that  an  axial  variation  in  the  focused  region  is  being 
observed.  These  variations  were  not  so  evident  in  the  plots  for  the  5.06  cm  and  3.3 
cm  focal  length  lenses. 

C.  MEASUREMET^  OF  BREAKDCWN  EJffSSION  FROM  ARGON 

Guided  by  theoretical  calcxilations  which  predicted  short  excited  state 
lifetimes  during  laser  Irradiation,  woiic  was  done  to  examine  the  'imlssion  from 
the  ionised  region  of  the  gas.  In  this  woric  a  nonochrooator  was  utilized  to 
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separate  the  spectral  wavelengths  of  the  emission  with  the  set-up  shown 
in  Figure  Careful  examination  of  the  rate  of  emission  during  the 

laser  irradiation  would  give  an  Indication  of  the  nvechanism  by  which  the 
ionization  grows  provided  that  a  sufficiently  sensitive  photomultiplier  was 
available.  An  analysis  of  the  emission  from  argon  ionization  region  revealed 
that  the  light  emitted  during  the  laser  irradiation  was  identical  to  that 
shown  in  Figure  ni-5,  and  was  independent  of  the  wavelength  within  the 
visible  spectrum.  Both  the  risetime  and  the  amplitude  of  the  emission  were 
dependent  on  the  amount  of  energy  into  the  gas.  Tne  risetime  of  the  observing 
system  was  limited  to  about  6  nanoseconds  by  the  Tektronix  585  oscilloscope 
used  to  record  the  photomultiplier  signal.  Measured  rlsetlmes  for  the  initial 
light  emission  were  as  fast  as  20  nanoseconds  for  high  energy  shots  when  the 
12.7  cm  lens  was  used  to  produce  breakdown.  This  risetime  increased  to  80 
nanoseconds  or  more  as  the  input  energy  dropped  to  near  threshold,  ’..’hen 
shorter  focal  length  lenses  (the  5.08  cm  ai;1  3.3  cm)  were  used  to  produce 
ionization,  somewhat  longer  risetimes  were  noted.  In  al  T  of  this  work  the 
pressure  of  the  argon  gas  was  700  torr.  At  this  pressure  no  signs  of  a  line 
spectrum  were  noticed.  This  is  consistent  with  the  results  expected  on  the 
basis  of  line  broadening.  Essentially  the  Ar  spectrum  is  so  dense  that  line 
broadening  of  the  type  expected  from  the  absorption  observed  would  produce  a 
continuinn.  The  same  measurements  nade  in  He  (to  be  discussed  below)  showed 
line  broadening  of  about  100  f.  during  the  time  of  laser  irradiation  of  the 
gas.  This  same  line  broadening  was  not  present  after  the  laser  irradiation 
turned  off. 
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Figure  III-4  Set-up  to  Measure  Spectral  Einission  ’/avelengthr> 
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Additional  Investigations  >fere  made  of  the  tljne  history  of  emission 
from  the  Ionized  volume  after  the  laser  radiation  subsided.  This  emission 
had  a  time  history  very  dependant  upon  the  wavelength  being  monitored.  The 

o 

emission  time  varied  from  about  15  microseconds  for  wavelengths  around  3^00A 

0  o 

to  over  200  microseconds  for  wavelengths  around  4600A  to  5600A  and  then  dropping 
off  to  about  20  microseconds  aroiind  630oX.  A  typical  time  resolved  photograph 
of  the  emitted  light  la  shown  In  Figure  III-6. 

D.  rOIIIZATIOi:  TIME  AS  A  FUNCTION  OF  FXAL  SPOT  DIA}!ETEn 

To  obtain  an  Indication  of  the  effects  of  the  focal  sp>ot  diameter  a 
series  of  measurements  were  made  using  the  cell  filled  with  700  torr  of  argon. 
Each  of  the  lenses  whose  focused  spot  dlamater  had  been  measured  previously 
were  used  to  produce  Ionization.  The  curve  shown  in  Figure  III-7  illustrates 
the  effect  of  spot  diameter  on  the  breakdown  (attenuation)  time  versus  peak 
electric  field  strength.  Repeated  experiments  have  shown  conclusively  that 
the  threshold  for  the  3.3  cm  lens  is,  indeed,  higher  than  the  thresholds 
seen  for  the  12.7  cm  and  5.08  cm  lenses.  It  may  be  that  the  higher  threshold 
Is  caused  by  an  Increase  in  the  diffusion  losses  which  accorapnay  the  smaller 
focal  di.ameter.  This  is  supported  by  the  fact  that  the  average  data  for  the 
5.08  cffi  lens  lies  above  the  data  for  the  12,7  cm  lens. 

It  should  bo  explained  that  in  the  course  of  experimentation,  it  became 
necessary  to  Increase  the  triggering  sensitivity.  This  caused  the  addition  of 
a  small,  constant  Increment  In  breakdown  time  measurements.  In  order  t'^  avoid 
confjslor.  all  of  the  data  presented  has  been  adjusted  to  produce  the  pulse  peak 
at  125  nanoseconds  from  the  start  of  the  pulse.  By  doing  this  data  pre¬ 
sented  In  earlier  reports  can  be  accurately  compared. 
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Fi/mre  III-^  >eakdov/n  ’^asion  fron  the  riasra  Fomed  \rron 
(20  u  s/cr.  horizontal;  2v/cn  vertical,  U500  ■) 
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Peak  electric  field  strength  used  for  Figure  III-7  has  been  calculated 


using  F 

P 

where 


(Fq.  ni-1) 


F  ■  peak  electric  field 
P 

P  input  power  determined  by  dividing  the  half-power 

pulse  width  into  the  total  energy 

Zo  ~  the  impedance  of  free  space 

and  A  the  area  of  the  half  intensity  contour  computed  using 

A  7d^  where  d  is  the  half-intensity  (3db)  diameter 
k  measured  for  each  lens. 


IE.\3lTRFIFUTrj  OF  ’"ni\KDO'./M  DI  HFLIUI^ 


A  glass  cell  similar  to  those  used  for  argon  was  constructed  and  filled 
with  700  torr  of  helium.  Attached  to  one  side  of  this  vessel  vms  a  small 
enclosed  capsule  of  mercury  sealed  so  that  it  could  be  opened  at  the  desir  ;d 
tune  by  a  small  ferrous  breaker.  Initially  the  mercury  capsule  i^mained  sealed 
and  measurements  of  breakdown  in  helium  were  made.  The  most  obvious  difference 
between  argon  and  helium  was  that  heliiom  showed  little  or  no  sign  of  attenuation 
of  the  laser  light  although  a  spark  was  visible  and  omission  was  measured  by 
the  monochromator  system.  Pecause  of  this  lack  of  attenuation  no  accurate  mea¬ 
surement  of  variation  in  the  time  of  breakdown  occurence  could  be  made. 

The  time  of  the  start  of  the  initial  light  emission  (from  start  of  laser  pulse) 
was  measured  for  a  largo  number  of  cases,  and  breakdown  as  measured  by  time  of 
light  emissioi^  seemed  to  occur  within  a  few  nanoseconds  of  the  peak  of  the  laser 
pulse  no  matter  how  great  the  change  of  energy  into  the  gas. 


’.Jith  available  energy,  breakdown  could  not  be  produced  with  the  12.7  cm 
lens,  .bo  breakdown  threshold  was  determined  by  using  the  5 .OF  cm  lens  and  the 
3.3cn  lens.  Knowledge  of  the  spot  size  of  each  of  those  lenses  enabled  the 
field  strength  to  be  computed. 


The  threshold  determined  for  helium  using  equation  III-l  was  7  x  10  volts/c- 
as  compared  to  3.9  x  10^  voltsAn  for  argon.  The  pressure  in  both  cases  vns 
700  torr;  however,  the  12.7  cm  lens  ^vas  used  to  determine  the  Ar  threshold  ;r.ile 
the  5.08  cm  lens  was  used  to  determine  the  He  threshold. 

An  attempt  was  made  to  determine  the  effects  upon  breakdo\m  caused  !y  a 

large  amo\mt  of  mercury  vapor.  To  do  this  the  mercury  ampule  attached  to 

the  cell  was  broken  and  the  mercury  vapor  was  allowed  to  r±x  wic-h  the  hcli’m- 

gas.  At  room  tenperature  and  below  no  effects  in  the  helium  breakdo^'m  were  o!  - 

served,  \fter  heating  the  ampule  to  drive  mercury  vapor  into  the  cell  it  ins 

possible  to  detect  a  noticeable  effect  on  the  breakdown  emission.  Ahereas  a 

line  spectrum  was  observed  in  pure  helium,  a  continuum  became  evident  when  r.crcur- 

was  added.  It  was  observed  that  when  the  field  strength  was  below  the  thresh. old 

of  helium  a  continuum  with  superimposed  Hg  lines  was  predominant  whereas,  v/her. 

the  field  strength  was  above  the  threshold  for  helium,  the  helium  lines  (notably 
0 

//^5A)  were  superimposed  on  the  mercury  continuum.  Any  tendency  of  the  nercur'^ 
ir’.p  -ity  to  lower  the  threshold  could  not  be  detected  above  the  emission  from, 
thf  Hg.  It  is  difficult  to  determine  whether  the  initiation  of  breakdovm  tak-'s 
place  in  the  main  gas  (Ho)  or  in  the  impxirities  present  in  the  main  gas.  The 
addition  of  Hg  to  easily  observable  gases  such  as  argon  may  give  more  indication 
of  the  (effects  of  impurities. 

Photographs  were  taken  to  record  the  time  history  of  the  breakdown  emission 
of  pure  helium  during  the  tine  of  irradiation  from  the  laser.  These  curves 
(an  example  shown  in  Figure  III-8)  compare  rather  well  with  those  from  the  break¬ 
down  in  argon  at  the  sane  pressure.  A  record  of  the  entire  emission  from  the 
breakdown  region  of  the  helium  is  shown  in  Figure  111-9.  No  explanation  can  be 
offered  for  the  shape  of  the  emission  c\irve  at  late  tines. 
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Portion  of  curve  shown  in  111-8  cannot  he  seen  in  this  cur/e 
because  of  its  sharp  rise,  short  duration,  and  high  amplitude. 
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F.  'lOME  0’’'''Fll7\'"TONr>  AT  LOJ  PRZZ3mr. 


\5  pressure  Is  lowered  the  diffusion  loss  from  the  breakdown  region 
omes  greater,  and  at  some  pressure  becomes  great  enough  to  prevent  the 
avalanche  buildup.  If  one  carefully  measures  the  number  of  charges  produced 
as  the  pressure  is  varied  below  the  avalanche  region  an  indication  of  the 
ionization  build-up  process  may  be  obt.vined.  In  order  to  make  measurements 

_7 

at  low  pressures  a  glass  call  was  constructed  and  evacuated  to  10  torr. 

This  coll  (residual  gas  air)  had  a  mercury  anqjule  and  the  planned  work 
was  to  increase  uhe  vapor  puressure  in  the  cell  by  increasing  the  temperature 
of  the  mercury  ampule.  Iloctrodea  were  provided  to  collect  the  charge 
produced. 

efore  mercury  was  allowed  into  the  cell  several  laser  shots  produced 
10^*^  charge  puirs.  After  several  shots  this  level  dropjjed  two  orders  of 

Q 

magnitude  ai*d  remained  at  10  charges  throxjghout  a  series  of  nhots.  Thereafter 

no  change  was  noted  when  the  mercury  was  allowed  to  enter  the  coll  and  raised 

_2 

in  vapor  pressure  to  10  torr.  The  charge  collection  current  was  different 
at  all  of  those  pressures  than  it  was  at  700  ton*.  The  current  had  three  peaks 
spaced  about  lOus  apart  and  increasing  in  amplitude  with  time.  It  is  not  yet 
clear  what  causes  this  shape  or  why  increasing  the  press^ire  of  Hg  had  no  effect. 
Perhaps  roflo''tod  light  from  the  rear  of  the  cell  produued  photo  electrons  from 
the  electrodes  which  were  collected.  This  work  Indicates  that  low  pressure 
investigations  may  produce  interesting  resxlts,  and  may  perhaps  provide  some 
answers  to  the  initiating  and  growth  mechaxd.8]n8  provided  more  care  in  experimen¬ 
tation  is  taken.  Low  pressure  work  will  be  continued  for  various  other  gases. 
Additional  precautions  to  eliminate  unwanted  variables  will  be  taken. 


G.  EFFECTS  OF  HIGH  AND  LO'J  PREGSUiIE 

A  test  cell  (See  figure  III-IO)  capable  of  withstanding  pressures  of 
1500  psl  tma  been  constructed  and  is  being  used  to  investigate  gas  ionization 
as  a  function  of  pressure.  Constructed  of  stainless  steel,  the  cell  is  capable 
of  being  evacuated  to  10~^’  nm  Hg  before  being  filled  with  the  desired  gas. 

A  Heise  guage  la  being  used  to  monitor  the  pressure  accurately  (jp.5  psi)  fron 
atmospheric  to  1000  pal.  An  additional  glass  cell  system  capable  of  evacuation 

-9 

to  10  nn  Hg  has  been  constructed  for  the  low  pressure  work.  This  cell  and 
associated  equipment  are  shown  in  Figure  III-U. 

H.  CONCLUSIONS 

The  process  of  gas  ionization  shows  a  threshold  dependence  upon  the 
intensity  distribution  in  the  laser  beam.  Since  significant  changes  in  the 
intensity  distribution  can  cause  significant  changes  in  spot  size,  care  must 
be  taken  to  monitor  these  variations.  Secondly,  spot  size  for  short  focal 
length  lenses  cannot  be  detemined  from  a  knowledge  of  beam  divergence,  thus 
flLCCiirate  measurements  of  spot  size  is  an  in^jortant  parameter  to  be  defined 
experimentally.  Reducing  the  spot  size  by  use  of  a  5.08  cm  focal  length  lens 
lowered  the  time  of  breakdown  for  a  given  field  strength  as  compared  to  that 
for  a  3.3  cm  lens. 

Results  of  the  work  done  in  helium  indicate  that  the  spectral  emission 
of  the  plasma  differs  for  the  various  gases.  The  easily  lonizable  mercury 
vapor  had  no  detectable  effect  on  the  breakdown  time  in  helium. 

I.  FURTHER  VroRK 

It  is  planned  to  investigate  the  entire  pressure  range  agreed  upon  using 
several  pure  gases  (Argon,  Nitrogen,  Relixia  and  air).  Work  is  also  anticipated 
to  determine  the  effects  of  a  neodymium  frequency  (1.06u)  pulse  on  several  gases. 
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A.  ni  'Rrr'CTI  jN 

1  2 

In  previ^  :3  iisc’issions  *  of  tho  '  reakriown  ~>f  f^ases  }y  la'-.ers  we 
have  conaiiiered  vari  I'is  nsnects  of  the  h  ■'pot^’.esia  that  at  hii’'  gae  nrcso'ires 
( one  half  at-oaphoro  ami  atove )  the  breakdown  is  the  result  of  initiation 
by  a  multi-plioton  ionizat'  in  event  follo;;ed  oy  an  avalanche  type  j^rowth  of 
the  ioniaation  to  tiio  observed  value.  In  this  section  we  will  consider 
further  the  avalanche  type  growth  of  tho  ionization.  We  will  aasune  that 
for  the  gaaea  under  consideration  an  ele.-tronically  excited  atom  is  rapidly 
ionized  by  multi-photon  absorption  so  that  tho  effective  rate  of  ionization 
is  equal  to  the  rates  of  electronic  excitation  plus  ionization  by  electron 
impact.  This  assumption  needs  to  be  Inveatigated  further  but  appears  to 
be  reasonable. 

Tho  topics  to  be  considered  in  this  section  are; 

B,  Further  evidence  conrerninf^  the  applicability  of  the  extension  of 
microwave  theory  to  tho  laser  breakdown  proi  lem  as  obtained  from  steady  state 
and  time  denen'ient  treatments  of  the  photon  absorotion  nrocess  using 
simplified  models  for  the  electron-gas  atom  coll isi ons, 

C,  The  calculation  of  avnlancfie  growth  constants  for  Ar  and  using 
the  extended  microwave  t'  eir;/,  .ond 

D,  Tiie  calculation  of  apnroximato  cic’ves  of  tl  9  thrt.sho,.!  laser 
intensity  required  for  breakdown  as  a  function  of  gas  density,  laser  pulse 
length  and  focal  spot  size. 


Thia  discussion  is  by  A.  V.  Phelps, 
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B,  iic  iivalence  of  Dijjcrete  ’’hoton  and  !!inrowavo  Theories 

Published  discissions  ^  of  laser  induced  breakdown  differ  considerai ly 
as  to  the  applicability  of  extensions  of  microwave  theory  to  optical 
frequencies  for  the  calc  ilation  of  the  rate  of  f^rowth  of  the  ionization, 
de  :ave  s-  own  oreviously  ^  that  the  differential  formulation  of  the 
Boltzmann  equation  used  in  the  n  crowave  calculations  is  a  pood  annroy.-a- 
tion  to  the  difference  form  of  the  Boltzmann  equation  anpropriate  to  the 
absorption  of  photons  in  free-free  transitions^ orovided  that  the  electron 
energies  of  importance  are  large  compared  to  the  photon  onergo^  and  that 
the  electron  energy  distribution  function  does  not  change  too  rapidly  in 
an  energy  Interval  equal  to  the  photon  energy.  In  practice,  this  means 
that  we  expect  the  microwave  theory  to  "ivo  good  values  for  the  rates  of 
excitation  but  possibly  to  yield  ionization  rates  which  are  too  low.  As 
further  evidence  for  the  correctness  of  this  prediction  we  shall  compare 
directly  solutions  of  the  difference  and  differential  forms  of  the  Boltzmann 
equation  using  approximate  forms  for  the  croaa  section  for  photon  absorp¬ 
tion  ancl  electron  excitat  ion,  dsing  one  s  ich  model  we  shall  inquire  as  to 
what  a  steady  state  solution  tells  us  al  o  t  the  rate  of  excitation  due  to 
the  at  sorption  of  a  finite  nunt  or  of  pnotons  as  comoar-  •  ‘o  the  r.ate 
calculated  using  ttie  microwave  ttieory.  This  investigation  is  largely  in 
response  to  ttie  suggestion  that  beca  ise  the  ionization  or  excitation 
potential  is  a  relatively  snail  multiple  of  the  photon  enerj",’,  the  electrons 
may  gain  sufficient  energy  to  excite  or  ionize  without  siffurini:  the 
randomizing  effects  of  large  numbers  of  collisions  implied  by  the  differential 
f  orm  of  the  Boltzmann  equation.  This  analysis  has  been  carried  out  using, 


random  walk  theory  hy  Drs.  .  Gaver  and  J,  A.  Marshal  of  the  I^aboratories 
'lathematics  Department.^  The  second  model  allows  a  solution  of  the  time 
deoendent  difference  equation. 

(A)  Steady  State  Solution 

The  lioltzmann  equation  on  which  the  statistical  analysis  is  based  is 

^q.  (6)  of  Ref.  1,  If  wo  assume  that  the  electron  energy  distribution 

6 

1  unction  consists  of  a  series  of  delta  functions  separated  by  the  ohoton 
energy  hr".  tlien  Sq,  (6)  of  Ref.  1  becomes 


/ pk4>^n 


where  the  index  i  is  an  integer,  is  the  fraction  of  the  electrons  with 
an  energy  of  ■  ih*)^  ,  and  are  the  probabilities  for  the 


absorption  and  stimulated  emission  of  a  photon  by  an  electron  vdth  an 
energy  ■  ihV.  From  Eq.  (5a)  of  Ref.  1  an  approximation  to  valid 

for  and  low 


Jl.  4^  ^ 
3  /if 


(2) 


2  2  2  -1 

where  £  •  e  E  (mw  )  is  the  classic*!  oscillatory  energy  of  a  free 

0 

electron  and 'T'  {(.,)  is  the  frequency  of  momentum  transfer  scattering 
'  m  •'i 

collisions  of  electrons  by  the  gas  molecules.  At  present  we  do  not  have 
acc  urate  expressions  for  and  so  will  consider  a  second  aporoximation 
to  which  is  probably  better  when  h  and  ar«  comparable.  This 

7 

relation  has  been  suggested  by  Holstein  and  is 
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(3) 


Li.  /£ .  ^  hJ^ ) 

^  h'y  X  ) 


In  order  for  the  electron  distribution  to  be  limited  to  the  discrete  set 

of  P.  's  the  energy  los.'-  in  excitation  f  is  ass-unod  to  be  equal  to  an 
i  ^ 

intef^ral  multinle  of  V*  i-'-  •  H.  From  the  detailed  balascinr; 

argument  given  in  Lq.  (7)  of  Ref-.  1,  we  have 


In  order  to  simplify  the  problem  we  will  assume  that  excitation  oci-urs 
as  soon  as  the  electron  reaches  the  level  H,  This  is  equlvalert  to  the 
assumption  of  an  infinite  excitation  cross  section  discussed  by  Allis? 

In  the  microwave  case  this  approximation  has  been  found  to  give  reasonably 
accurate  excitation  rates  at  electric  field  intensities  which  are  not  too 


The  frequency  of  exciting  collisions  using  the  model  described  above 
is  equal  to  the  reciprocal  of  the  expectation  value  of  the  first  passage 
time  \  This  quantity,  E(TIH),  is  given  by  ^ 


Equation  (5)  has  been  eval’uated  as  a  function  of  H  for  various  assumed 
dependence  of  In  this  discussion  we  will  consider  three  cases;  ^ 
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a) 


^  hy  /ly 


(6a) 


Case  a)  c^n■esn''.d3  to  the  sul  stitution  in  Eq.  ( 3 )  an  energy  Independent 

freqiency  Tf  nonentum  transfer  collisions  such  as  found  for  electron- 

neutral  collisions  in  and  M  at  enorpies  above  a  few  electron  volts. 

2  e 

This  is  expected  to  he  the  rmst  realistic  of  the  models  conside-red  here 
since  electron-neutral  collisions  dominate  the  nhoton  ahsorotion  process 
d  irinp  the  fjro.Jth  of  ionization.  '.<o  note  that  the  use  of  the  constant 
collision  fr'-'juency  assumption  in  Eq,  (2)  for  our  simple  model  leads  to 
■  0  and  trapping  of  electrons  v.l  zero  energy  and  so  is  unrealistic. 

Cases  b)  and  c)  correspond  to  the  assumption  of  electron  scattering’  by 
ions  in  which  the  collision  frequency  is  approximately  inversely  propor¬ 
tional  to  the  3/2  power  of  the  electron  energy,  Case  (c)  is  of  interest 

because  it  leads  to  A  ■  QO  and  because  it  is  a  case  which  one  w  )ij.d 

o 

expect  to  find  treated  by  others.  In  addition  for  case  (c)  X  ^  -  JC 

in  Eq.  (1)  so  that  the  electron  motion  in  energy  is  a  pure  random  wal'’. 

Rather  than  show  the  calculated  excitation  rates  directly,  we  will  compare 
them  with  the  results  of  the  differential  or  microwave  t‘  eory. 

The  sol  ition  to  the  differential  form  of  the  Boltzmann  equation  for 
the  single  m  idel  considered  here  is  obtained  by  setting  the  right  hand 
side  of  Eq,  (9)  of  Ref.  1  to  zero,  solving  for  fi^)  and  then  normalizin;’:. 


i 


IL  can  lo  shown  that  this  procoriurc  leads  to  the  relation 


%  ~  i(hy) 


j  clC 

MO 


(7) 


for  case  (a),  XCc)  =  z  ) 


and 


£ 

For  case  1),  \(l)  —  ('£•  +  ) 


and 


^  -  ^\ihy  j 


For  case  (c)  \(i)  =  At  C^y'/s) 


and 


-  Jt_  /  O  N  ^ 

fy,  If  hjlJ 


(  ''•■'I  ) 


(9a) 


(%) 


(9b) 


(Ic) 


(9c) 


The  ratio  of  the  excitation  frequency  calculated  using  Eq,  (5)  to 

that  given  by  Eqs,  (9a,  b,  c)  is  s  own  in  Figure  IV-1,  The  lower  curve 

shows  that  for  the  model  of  electron-neutral  scattering  in  which  the 

12 

electron  collision  frequency  is  indeoendent  of  energy,  the  excitation 
freq  lency  using  the  discrete  nodel  is  lower  than  that  calculated  for  the 
continuous  model  by  roughly  In  the  model  of  electron-ion  scattering 

the  dashed  curves  show  that  the  calculated  excitation  freq'iency  for  the 
fiiscrete  case  is  generally,  but  not  always,  higher  than  that  calculated 
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Ratio  of  Excitation  Frequencies  —  (discrete)/ (  con' t) 


Ly  the  c  )nLin  iTis  case.  I'or  _>  the  diffei-- T.ce  i;;  lens  t:.an  ,i . 

Jince  these  calcalati  :)n3  sh'jw  that  &hc  crnr  resultiny  fnn  the  .se  )!’  t-.e 
cor;tin'Jou3  theory  ilopunda  stron^^ly  on  the  er.eryy  deyendeiice  of  the  al  s  r  h  or. 
irobability,  it  in  very  desirable  to  carry  out  calculations  of  t.ue  rates 
of  excitation  usinp  nore  realistic  cross  sections  usiru:  the  liifferer.ce 
forr.  of  the  lloltznann  equation,  such  a  caiculation  is  bein,;  atteT.p'  ed  \  ut 
so  far  has  leen  unauccessf ul , 

(b)  Tire  Depenaent  solutions 

As  an  additional  check  on  the  applicability  of  the  extended  nicro- 
wave  theory  to  the  prediction  of  the  rate  of  growth  of  ionizati''n  due  to 
the  absorption  of  ohotons,  v'o  liavo  calculated  sonic  tine  ieoen  'er.t  el<!ctr"'n 
e-.ercy  distribution  functions  and  excitation  f r<'q  icncics.  These  -"I  es 
can  then  be  compared  vrlth  the  results  of  the  usual  microwave  the 'ry  ich 
assumes  that  electron  energy  distribution  is  indeoendent  of  tine.  7hc  ".o  ;ol 
used  in  these  calculations  is  more  exact  than  that  discussed  above  in  that 
t..e  probability  of  an  excitation  collision  is  finite  ,  altuou -h  of  s  nple 

13 

analytical  form.  As  above,  we  neglect  the  energy  loss  in  elastic 
collisions  with  the  gas  atoms,  although  t.his  loss  is  ex])ected  to  be 
important  at  high  pressure  in  the  rare  gases.  ^  The  ocimtions  which 
have  been  solved  are: 

H  oo) 


ill 

dl 


I  7 


whore  i  -  V  for  i  2-  ^  -  0  I'lr  i  ^  Ii;  -  yC  /L  Tor 

^  M  j  j  «r  (^-0  z  /  for  X  z  !•  -^=0  for  «5// 

id  ~  ^  ^  r^'4-  N’ote  Lhnt  ;iere  we  luivo  u-^otl  the  f  irn  if  the  con.;tant 


c  illiiiioii  frequency  a  roxirvitioii  correapon  i^n  to  r,q,  (2).  Thia  ainpler 

r  ^ 

approxiraution  '  ia  jioaaible  aince  with  llio  preaont  noilol  inelaatic  col- 
liaions  ana  atinulatod  photon  onisaion  do  not  oori'olato  the  level  at 

“  0.  ..'o  note  that  the  ays  ten  of  eq  iatlona  is  terr-'i.nated  .t  i  =  !!, 

where  i;  is  chosen  by  trial  anci  error  1 1  be  large  enough  so  that  the 
nopulations  at  energies  of  interest  are  ind.ojiendent  of  M,  >l6  can  use 
the  calculated  val  ics  to  deternino  the  excitation  freq  ency  Kiinf;  the 


relation 


=  2. 

i--n 


rinally  we  note  ttiat  o  ir  asa  inption  as  to  the  energy  denenderico  of  li. 


is  eq  .ivalent  to  assuming  that  the  cross  section  lor  excitation,  is 
given  by 

/Vv  = 

7  1/2 

wtiere  v  ■  5.93  x  10  (£)  cn/sec.  For  He  and  electron  energies  within 

a  few  electron  volts  of  the  excitation  threshold  H  ■  11,  and 

^  l.ii  X  10”^^  sec”^.  The  calculations  given  below  have  been 
carried  out  by  D.  Wei  of  the  l^aboratories  Comouter  Oro  ;p. 

The  results  of  the  solution  of  Eqe.  (10)  are  s'rwn  in  Figure  .'  .-2 
where  is  plotted  as  a  function  of  i  -  £^/h')/for  B  ■  15.  These 


h 


S3 


rnsults  are  for  the  case  in  v;':ich  ■  1  at  t  a  0,  !!ere  tine  t  is  ncns'.ireci 

in  units  of  1/  sic',  that  T  t.  >!e  sec  that  the  Histrii  .tion  furiCti  'n 

has  r' ac.ued  its  steady  val  :e  i  y  T  ■  6.  'Jslnc  iiQ.  (2)  and'Ty  ■  7.2  x  10 

cn^/sc“c  as  expected  for  helium  at  a  g&s  dcr.s.  ty  I.',  the  c'lrves  in  Fijp.iri' 

rV-2  corrns'iond  to  the  case  of  a  field  stro:,[;th  d  ■  3.’'  x  10^  v/cn  at  the 

o  iJ 

laser  freqiericy  corresp  ndiiip  to  A  or  a  laser  poi;er  of  j,  I  x  iJ 

o  19 

watt/cn"  a:;il  to  a  pas  derui .  *.y  :>f  2.1,5  X  10  atons/cc  or  ono  atmosphere 

at  30  7°K.  In  this  case  A  ■  2.3  x  10^  soc  so  that  the  distribution 

-9 

function  has  reached  its  final  form  by  about  2.5  x  10  sec. 

In  the  present  I'roi  len  we  are  primarily  concerned  witli  the  time 

depeiulence  of  the  excitation  frequency  as  shown  In  Fif^ire  IV-3.  The  lov/er 

two  curves  (P  •  l5)  s!  ow  the*^/^  values  calculated  usinf^  iq.  (11)  for  the 

model  discussed  above  for  helium  and  for  two  different  assumptions  as  to 

the  initial  state  of  the  electron,  i-e.,  p  ■  1  at  t  •  0  and  P  •  1  at 

6  i 

t  ■  0.  i.'e  see  th't  the  excitation  frequency  has  reached  its  final  value 

-R 

of  2.6  x  10  sec  Iv  about  2.5  nsec.  The  upner  curve  shows  the  excita¬ 
tion  frequency  calculated  for  b  ■  0.5  cor  osponding  to  helium  at  r‘nos- 

7  ,  12 

p'.oric  press  ire  and  to  a  laser  field  strength  of  2  x  10  v/cm  or  l.O^i  x  10 
2 

watt/cm  .  T  is  calculation  has  not  been  carried  to  completion  Leca  :se  of 

the  largo  amount  of  compute-r  time  r-'quired  but  it  ap;)oars  that  the  evcita- 

-10 

tion  rate  at  10  sec  is  sit'nif icantly  below  the  final  value.  The  time 
required  to  attain  the  steady  stave  excitation  rate  is  important  I'ecause 
each  new  electron  in  the  avalanche  will  have  to  indorgo  such  a  transient 
period.  The  average  excitation  rate  is  thei*efor''  an  ap-'ropriato  integral 
over  the  instantaneous  excitation  rate  weighted  according  to  the  pro¬ 
bability  of  finding  a  new  electron  in  a  given  initial  state  after  the 

photo .^nization  ^f  the  excited  atom.  Je  have  not  car-'ied  uit  this 
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Electron  Energy  in  Units  of  hi/ 


Figure  IV  -  3  Tine  Dependent  Ellectron  Energy 
Dietributlons 
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calc  1  fit.  ion. 


Ttiir.  avora/’o  ''xc  ita^  ion  frnqicnry  vrlll  fv-iinrally  be 
sriallnr  t'.an  ttie  final  val  in  beca  u;g  ti-  s  calculation  is  one  way  of 
allowin/;  for  tbe  onerf^^  rnq  iirn‘H  to  r.oi  ie  ‘.he  nnorpy  of  each  now  electron 
to  '  tio  mean  electron  en-nty. 

.j’e  now  wish  to  compare  ttie  final  excitation  rates  with  those  calc  ilated 
isinr  oth'  T  aporoximate  midels.  Isind  ttie  steady  state  solution  corres¬ 
ponding  to  the  tine  dependent  calc  .lations,  i.o.  dq.  (2)  with^  ind.e-.endc.it 

^  n 

of  on  rrj  and  the  summations  in  nq.  ([>)  i;oin;',  from  1  to  •!,  we  find 

,  .  -1 
■  1.1  X  1  )  conpaiMd  to  a  fxrvol  value  of  1.2  ^  0.‘,>  x  1-  usiny  tin; 

tine  dependent  silution  for  I  =  l^i.  The  agroement  botwoon  ttiose  results 

is  expected  since  they  are  based  on  models  wtiich  are  neyilipibly  diffeent. 

7he  Comstar. t  collision  frequency  model  usinp;  Kq.  (3)  in  dq.  ('^)  yields 

fx/X  •  I.IJ.  x  10"^  d.ue  to  the  slightly  smaller  values  of  ?-i/X  over 
ttie  ran,:e  i:  i  considered  )  to  !I  rather  that  1  to  If  +  1),  The  soluti  on 

to  tfie  difl'erential  form  of  the  Poltzmann  equation  cives"^^/  ^  ■  1,^7  x 

-1  -2 

10  and  1,5'^  X  10  w  en  dqs,  ''2)  and  (3),  resnectively,  are  sulstitutori 

into  dq,  ( 2 )  and  X  is  constant.  The  extended  microwave  tfieo/y 

1  3 

using  realistic  cross  sections  as  tiiscussed  in  iief,  1  gives"^  ■  1,5  x  10 

X  ^ 

sec  ”  or  "y/  “X  ■  ^  10  •  This  lower  value  is  largely  the  result 

'  X 

of  the  neglect  of  energy  loses  due  to  '^lastic  collisions  in  the  simple 
models  we  have  considered  in  this  section.  Finally,  we  note  that  tiie 
steady  state  value  of  for  B  ■  0,5  is  greater  than  9  x  10  ^ 

-  6.1  X  10^  sec"^  or^^  X  “ 


as  compared  to  the  microwave  values  of 

-2 


9,7  X  10  ,  In  this  case  elastic  collisions  are  much  less  important  than 

in  the  lower  electric  field  (higher  P)  case.  These  calculations  again 
suggest  the  possibility  of  significant  differences  between  steady  state 


h.- 


excitaLiTri  rates  jsin-  tho  twi  I'or'.j  of  the  i'olt  "oan  fcuation.  r, 

i  cy  S'.ow  that  a  r  ally  valic  c  oroarison  cr-n  )o  rv  de  only  let*"!' 

n  )di;l.;  ar*'  lood  to  ac -ount  ir  t-.e  effects  s  c:;  as  olavO^c  c^lli.';.  ir.r. 
a:.'!  t’le  energy  r-'q  ;iro  1  to  heat  op  the  now  cli  ctr iniJ . 

J.  Avalanche  Irowth  in  Ar  and 

I:.  t,.'.3  secti  in  we  precent  toe  results  of  "-ale  ilati  onj  nf  cieff  ^  3 

ynvm.in,-  the  rate  if  pnwth  of  ionizati  ri  in  A.r  and  i.'p  ani  in  the  caoe  i; 
Ar  '..'e  apply  theje  reniilts  to  tht  prediction  ^f  the  tine  requir' c  for  Iri’ak- 
(iiwn  in  /or  'on  lor  sinplified  experincntal  conditions.  Tho  notati  .n  an  i 
sinp'lif ications  here  arf-  identical  with  t.h.at  of  .tef,  1  and  will  not  he 
re MO’  ted  here, 

Fi(Pire  rV-u  shows  values  of  ttio  excitation  and  ionization  free^  oncies 
per  atom  arifi  tho  electron  diff-ision  coefficient  at  unit  dei.oity  for  Ar  at 
free  .encies  such  thatC</^>^.  These  calc'olations  W'  re  made  usinr,  tiio 

17 

cross  secti-ons  and  computer  nroyran  discussed  }y  Kngelhardt  and  ihelns.  ' 
Firp;re  IV-5  sh  iws  the  calculated  excitatior.,  ionization,  and  diffusion 
freq  encies  for  electrons  in  Ar  at  one  atnosphore  on  ssure  at  300°?!  as  a 
function  of  the  field  strenfth  of  the  ac  field  at  the  rrpy  oncy  of  a  rdy 
laser  (69li3  X)»  These  curves  are  t’.en  us^d  to  calculate  the  curves  s  iwn 
in  Fiyure  IV-6  of  the  tine  requir-'d  f  or  ".he  ionization  to  prow  iy  an 
avalanche  process  from  1  electron  in  the  focal  voIuitb  (  ICc^  olectron/c"- ) 
to  complete  ionization.  As  discussed  in  iief.  1  the  soliii  cm'ves  are 
calculated  under  the  assumption  of  electron  loss  ly  free  diffusion  to  tr.o 
cylindrical  bountlary  of  the  f  ocal  voliane  while  the  ^'.ashoii  curves  c  rres  ond 
to  nepligible  diff  ision  loss  as  when  antipolar  diffusion  is  in;  ortant. 

;‘9 
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Figure  IV  -  5  Excitation,  Ionization  and  Diffusion 

Frequencies  for  Argon  at  One  Atmosphere, 
3OOOK  at  69U3  A  Using  Microwave  Theory 
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Fijjure  IV  -  6  Time  Required  for  Laser  Breakdown  In  Argon 

at  One  Atmoaphere  and  10”  *”00  Radius  Focal  Voijr 


It  ia  still  rathor  risky  to  con’virc  the  ’•  .-iults  of  t'.:;;  t‘ e  ry  '..'ith 

ox-u'r inent  lecaiso  of  the  uncertainties  in  the  tine  o!'  initiation  if  the 

avala;.c!;o,  the  snatial  and  time  variations  in  the  laser  nower  (ie-.r..t;.' 

w-  .c)i  are  ■  ot  included  in  the  theory,  ant'  the  uncertainties  remain  i  r. 

the  annlicalility  of  the  extentlo'!  nicrov;ave  toeory.  ' 'o;;ev  u’,  t.'e  r.  )te 

t  .at  v.’.e  threshold  fioli  strength  s'.o’.;:.  in  '.'i  ’lr':  17-6  f''r  freit  “lectrin 

difl'  ;sion  and  excitatini  follou'ed  ly  iihotoionizati  i.e.  }.2  •/.  10^ 

•..af/c  i''  or  3.’  x  10^  V/cn,  is  in  rather  good  agreement  ■..’.tti  ..’a;,nu'i:.t ' s 
10  2 

val  .‘i  of  Ii  X  10  watt/cn  for  a  SO'  large.'  focal  ratii  ;;;  I  ut  a  ..li  ntly 


lo'.;o;'  gas  lonsity. 

figure  IV-7  snows  calculated  e;-:citation  and  ionisation  fre'']"iCnc'’ es 

per  molecule  for  electrons  in  'onder  conditions  such  thritCt/ 

These  calculations  wore  made  using  the  cross  sectioixs  and  con  utationol 

19 

procedures  discussed  by  dngelhardt  and  Phelp.s.  ^  ,Jp  have  not  calcul:  teti 

the  excitation  frequency,  etc,  for  i;^  at  atnosnheric  pressure  since  no 

measurements  are  available  of  the  tine  of  breakdown  as  a  function  of 

laser  po\;er.  h'itliout  such  data  we  cannot  be  sure  of  the  sigrdficanco  of 

20 

the  observed  threshold,  A  preliminary  attempt  to  cornare  theory  and 
ex-'^erime  t  is  .’^Iven  in  th.e  next  section, 

I),  Approximate  Threshold  Predictions 

In  t!.is  section  we  will  make  of  the  exoitatio'  ,  ionization,  and 
diff  ision  coefficients  for  He  and  Up  which  have  been  given  in  i(ef,  1  and 
in  the  preceding  section  to  predict  approximate  rolati  ns  for  the  intensity 
required  to  produce  breakdown  of  these  gases  by  a  pulsed  ruby  laser,  ..'o 
shall  calculate  the  laser  intensity  necessary  to  satisfy  two  breakdown 


criteria.  I'irst,  we  assume  that  the  br  ar:down  threshol  d  is  dctfin  nod  by 

the  condition  that  the  electron  ionization  frecuency  ofi  ial  th<;  free 

difi'usion  frequency.  This  static  or  aiffusion  controlled  condition  is 

the  sane  as  the  conventional  microv/ave  br"akdo«n  criteria  and  has  been 

applied  l  y  several  authors  ^  to  the  laser  breakdown  nrol  ler..  Secondly, 

we  ass  tne  that  the  breakdown  thres'.ol  i  is  deternineri  by  t}ie  rocu  ronent 

t  at  the  ionization  Increase  from  one  electron  in  the  f  -ical  vol'ine  to 

co'iolebe  ionization  during  the  tine  the  Ipser  nulso  is  on.  'n  the  results 

shown  in  Figures  IV-c  and  IV-9,  we  have  continued  to  make  the  simplifying 

assumptions  discussed  in  Ref.  1,  e.g.,  a  spatially  uniform  laser  intensity, 

completely  absorbing  boundaries  at  the  focal  spot  radius,  a  constant 

amplitude  of  laser  intensity  for  the  duration  of  the  pulse,  negligible 

effect  due  to  electron  attachment,  and  the  presence  of  an  initiating 

electron  at  the  beginning  of  the  laser  pulse.  It  will  be  noted  that  we 

have  not  shown  a  limiting  curve  corresponding  to  the  threshold  laser 

intensity  required  to  produce  an  initiating  electron  during  the  laser 

2 

pulse.  As  indicated  previously  ,  one  can  predict  the  depeudence  of  this 

minimum  laser  intensity  only  if  one  knows  the  nature  and  density  of  the 

atoms  which  supply  the  initiating  electron.  Attempts  to  measure  this 

21 

ionization  rate  experimentally  and  to  nredict  it  theore* ically  have 

22 

teen  unsuccessful. 

The  static  or  diffusion  controlled  breakdown  criteria  is  obtained  by 
setting  It  •  0  in  Eq.  (13)  of  Ref,  1  so  that 


y,  -  y.  -- 


^*5 


(13) 


Figure  IV  -  8  I/iser  Breakdown  in  He 
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Figure  r/  -  9  I^ser  Breakdown  Thresholds  In  N2 
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Nitrogen  Density  ( mole cu les /cc  ) 


.  in^e 


wavo 


(V  +  V  )  /:;  and 
V  n^ry  unod  tiorc, 


'  !i  are  functions  or.ly  of  £/’I  and^yf!  in  the  micro- 
it  is  convenient  ti  rev/rito  !-q,  ''13)  in  the  f  im 


yvVl?- 


(11;  ) 


'lore  cor- esnojK^s  to  the  ns  .  1  ac  i nriizat i '>n  coef  1' ic'>'nt, 

hl-.en  UJ  neq'nirf  !  1' "yr  1  :■(•,•; h'io.-.n  c  n  !e  ict'nn.ineri 

2 

directlv  llron  a  plot  of  '7^  c’U’ves  fir 

i'in'U’ea  I.’-  and  17-9  s..nw  the  result  o;  s  ich  calculations  for  He  and  N2 

-2 

v;hon  li  is  taken  to  be  1C  cm.  The  fliflusion  controlled  thri.'Saold  curve 
Shawn  for  I;'.,  is  low  by  as  much  as  20'  at  tne  highest  der.jities  because  tne 

C 

corrections  due  to  a  finite^/i/  had  not  been  c.alculateri  at  the  tine  the 

23  ^ 

graph  was  pr  . pared. 

The  thr-'-sliold  intensity  required  to  satisfy  the  density  buildup 

2 

condition  is  calculated  ''ron  hq.  (13)  of  iief,  1  by  neglecting  the 

30 

term  and  ass'iming  that  the  electron  density  Increases  by  e  during  a 

60  nsec  constant  amplitude  laser  pulse.  This  relatively  large  value  was 

1  ’ 

ised  in  order  to  annroximate  conditions  in  the  exneriner.ts  of  '..'aynant 
rat- *’r  t!  an  those  of  fleyerand  and  H«  ;ght.  ^  Thus,  the  threshold 
condition  for  the  si.allower  curves  of  Figures  IV-d  and  h/-'^  is 


N 


H  ^  sec 


(15) 


Again  the  curve  shown  for  has  not  ywt  been  corrected  for  finite^/-^^  ^ 
and,  in  this  case,  is  low  by  as  much  as  a  factor  of  2  at  the  higher 


l*r 


iicMi3.  ^ies. 


./e  shall  concl'i'ie  this  discussion  i.ith  con-.cntc  ren'-rdinf^  tho  con-r 

ison  of  the  experimental  data  and  theoretical  curves  of  Fiipuojs  IV-  ■  and 

17-9.  First  v;g  note  that  thiC  two  threshold  c'lrves  cross  at  nos  donsit,i''s 

somewhat  below  atnosplieric  so  that  diffusion  losses  dotemine  the  *;■  res  ol  ; 

at  the  lower  gas  densities  and  the  buildup  of  electron  density  det''mirv!S 

the  threshold  at  the  highei  dor.sities.  Those  threshold  Return,  nin,-  c  u'ves 

are  s  lown  by  solid  lines.  In  the  case  of  V**,  .'i,;  re  IV-  '  shnws  t!/  •  t  .s 

1'’ 

rough  theory  is  in  reasonable  agreonont  with  the  data  of  ..aynant.  ..hen 
consideration  is  given  to  the  shorter  laser  p'ulse  l-nu.  ;th  used  ly  eyeran  : 
and  !!au  ^ht  the  agrttnent  is  reasonable  at  near  atnosphcric  Icns.ty  Vut  the 
t.'.eoretical  cur’/e  seems  to  decrease  too  slovjly  with  incr-;asi:.‘-  -as  dons. tic 
li.  It  is  to  be  noted  that  the  relative  values  of  the  br' akdov  thr''Shold 
in  de  given  by  I'inck  ^  arc  similar  to  those  of  I'eyor  nd  an'i  ’'au  •  t  in 
Figure  r/- 

Figure  IV-d  s'lows  a  comnarison  of  the  approximate  theoretical 
tl.rnshold  laser  intensities  for  H2  with  the  experimental  data  of  ’’eyerand 
and  Maught.^  This  comparison  is  reasonable  from  the  point  of  view  of  the 
microwave  t  eory  since  the  measured  Townsend  ionization  coefficient  for 

pi 

h'2  and  air  are  nearly  the  same  and  t’.eir  d-ffusion  coefficients  are 
exoected  to  le  similar.  As  in  the  case  of  helium  the  agreement  is  quite 
good  at  the  lower  gas  densities  when  one  considers  the  shorter  pulrc 
length  used  by  lieyeranl  and  Haught.  The  relatively  ra  id  decrease  in 
experimental  thros'uold  laser  intensity  with  increasing,  gas  de.-isity  c  omper*.' 1 
to  the  theoretical  curve  is  similar  to  tiiat  s.nown  in  Figure  IV-d  for  Ho, 
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Tliis  behavior  is  also  seen  in  the  relative  values  niven  by  I'inck, 

finiteO/y^  correction  referred  to  would  increase  the  discrepancy  even 

further  at  hip;h  f;as  densities.  In  both  He  and  it  would  appear  that 

the  decrease  in  the  efficiency  of  conversion  of  photon  energy  into 

excitation  (and  ionization)  predicted  by  the  microwave  theory  at  low  f/yis 

too  large.  In  the  He  case  this  energy  goes  into  elastic  recoil  collisions 

while  in  tne  Ng  case  this  energy  results  in  the  excitation  of  relatively 

19 

high  vibrational  states.  Obviously,  further  investigations  of  this 
behavior  are  necessary. 

E,  Sumiiary 

This  report  presents  evidence  that  the  absorption  of  relatively  high 
energy  photons  in  free-free  transitions  may  in  some  cases  lead  to  snail 
(  30, j)  but  possibly  significant  deviations  in  the  rate  of  growth  of 

an  electron  avalanche  from  that  predicted  by  an  appropriate  extension  of 
the  microwave  tlieory,  l.e,,  the  continuous  form  of  the  Boltzmann  equation. 
This  evidence  is  obtained  by  solving  the  difference  form  of  the  Poltzman” 
equation  in  both  the  steady  state  and  time  depe;xient  cases  for  various 
simple  models  of  the  electron-atom  or  electron-ion  collision  orocesses. 
However,  since  fliroction  and  magnitude  of  the  deviation  denend  upon  the 
model  used  and  since  several  important  affects  are  not  included  further 
comparisons  of  the  two  theories  are  desirable. 

Die  report  also  presents  additional  calculations  of  the  excitation, 
ionization,  and  diffusion  coefficients  and  makes  use  of  those  to  comuute 
avalanche  growth  times  and  approximate  ciu:ves  of  the  laser  intensity 
required  for  breakdown  as  a  function  of  gas  density.  In  general  the 
agreement  between  theory  and  experiment  is  within  a  factor  of  2  in  laser 
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intensity  or  UCXj  in  laser  electric  field  intensity.  There  are,  hov;cver, 
systematic  deviations  between  theory  and  experiment  which  should  be 
explored  further. 
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APPENDIX  A 


SOME  CONSEQUENCES  OF  LASER  BREAKDOWN  INITIATION 
BY  A  NON-LINEAR  PROCESS 


A.  V.  Phelps 

WesClnghouse  Research  Laboratories 
Pittsburgh  35,  Pennsylvania 


The  model  of  the  laser  breakdown  process  which  appears  to  best  fit 
the  presently  available  data  is  that  in  which  the  initiating  electron  is 
produced  by  non-linear  or  multiphoton  ionization  of  an  impurity  or  the  main 
gas  and  is  followed  by  an  exponential  growth  of  ionization  resulting  from 
either  electron  excitation  followed  by  photoionization  or  direct  ionization 
by  electron  impact.^  The  purpose  of  this  memorandum  is  to  examine  some  of  the 
consequences  of  the  assumption  that  a  non-linear  or  multiphoton  process  provides 
the  initial  electron. 

13  Time  dependence  of  ionization 

The  time  dependence  of  the  ionization  resulting  from  a  process  in 
which  the  photon  flux  acts  directly  on  the  neutral  atoms  is  significantly 
different  from  that  in  which  free  electrons  absorb  the  photon  energy  and 
produce  ionization.  The  continuity  equation  for  the  electron  density,  n,  in 
the  general  case  is 

"ahl  +  V-n-DV^n,  (1) 

t  xi  ' 

where  a  is  the  rate  of  ionization  of  neutral  atoms  or  molecules  of  density  M, 

is  either  the  frequency  of  excitation  and  ionization  by  elecr^-on  impact,^  and 
>)  is  the  electron  diffusion  coefficient.  Unless  otherwise  stated,  we  will  assume 
that  the  laser  power  increases  abruptly  to  some  constant  value  at  t  >  0,  so  the 
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a,  V  . ,  and  D  are  constants.  As  before,  we  approximate  the  solution  to  this 
XI  2  2 

equation  by  replacing  DV  n  by  Dn/Ji.  where  A  is  the  diffusion  length 
corresponding  to  the  fundamental  mode  in  a  cylinder  with  a  radius  equal  to  that 
of  the  focal  spot.  Equation  (1)  now  reduces  to 


.  a  H  ■.  -  v„)  n 


(2) 


where  ■  D/A-  .  The  solution  to  this  equation  for  the  case  of  n  ■  0  at 
t  ■  0  and  n  «  M  is 


n 


a  M 


II  -  exp  (v 


xi 


Vq)  t)  . 


(3) 


For  V  ,  <  v^,  this  solution  reaches  a  steady  state  value  of  ctM  (v_-v  ,) 

xi  D'  ^  D  xi'  ' 

whereas  for  it  increases  exponentially  with  a  growth  constant  of 

-  Vjj.  For  -  Vjj  ■  0,  n  ■  CKMt.  Note  that  our  assumption  of  n  «  M  may 
be  violated  at  rather  low  n,  if  M  is  an  impurity. 

The  mathematical  approach  which  is  used  above  does  not  show  the 

possible  Jlscontinuous  nature  of  the  true  breakdown  problem,  i.e.,  Eqs.  (1)  - 

(3)  are  valid  only  when  there  is  a  reasonably  large  number  (e.g.,  10  or  more) 

of  electrons  present.  This  point  is  important  in  evaluating  experimental  data 

since  it  means  that  if  the  rate  of  electron  production  by  the  non-linear  process 

is  small  enough  it  will  dotermine  the  minimum  laser  power  at  which  breakdown 

2 

is  normally  observed.  Thus,  if  the  effective  pulse  length  is  /.t  and  if 

19  5 

a  M  A  t  <  1,  then  even  though  (v^  ^  -  Vp)  At  >  In  (10  /lO  )  so  that  an  avalanche  could 

ionize  essentially  all  of  the  atoms,  breakdown  may  not  occur.  If  this  were  the 

case  breakdown  might  occasionally  be  observed  at  laser  powers  significantly 

below  that  for  which  a  M  At'^l,  i.e.,  breakdo%m  might  be  initiated  by  the  truly 

random  background  electron.  The  occurrence  of  a  single  rare  breakdo%m  event 

3 

at  power  levels  well  below  the  normal  threshold  has  been  observed. 

If  the  production  of  an  electron  by  a  non-linear  process  does 
effectively  determine  the  threshold,  then  observations  of  the  time  dependence 
of  the  growth  of  ionization  and  excitation  should  show  that  growth  constant 
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^  ^  4  ~  approach  zero  as  the  breakdown  threshold  is  approached 

but  that  the  observed  delay  is  due  to  a  delay  in  the  time  of  initiation  of  the 
avalanche.  Since  it  now  appears^  likely  that  measurements  of  the  light  intensity 
will  follow  the  ionization  build  up,  it  may  be  possible  to  extrapolate 
measurements  of  the  ionization  growth  back  to  the  time  of  appearance  of  the 
initiating  electron.  Such  measurements  will  require  the  use  of  photodetectors 
with  sufficient  sensitivity  and  speed  to  measure  the  light  generated  during  the 
portion  of  the  growth  period  when  the  degree  of  ionization  of  the  gas  is  small, 
e.g.,  when  the  light  output  is  ~1^  or  less  of  the  peak  value. 


2)  Dependence  of  ionization  on  geometry,  laser  power,  etc. 

The  requirement  that  the  laser  produce  the  initiating  electron  means 

that  in  general  two  distinct  threshold  conditions  must  be  met.  These  two 

conditions  will  have  rather  different  dependences  on  laser  power,  focal  spot 

size,  gas  density,  and  laser  frequency.  The  dependence  of  the  avalanche  growth 

condition  on  the  experimental  parameters  has  been  discussed  in  detail  for 

helium^  using  the  microwave  theory  as  a  basis  for  obtaining  quantitative 

predictions.  In  general,  this  analysis  shows  that  the  threshold  condition  for 

gas  densities  of  the  order  of  atmospheric  and  reasonably  small  focal  spot  radii 

is  determined  by  balancing  an  effective  ionization  frequency  which  varies 

approximately  as  N(E/u))^,n*>2  to  ^against  a  diffusion  frequency  which  varies 

2  -I 

approximately  as  (NR  )  .  Thus,  the  avalanche  threshold  condition  is  that 

(NR)^(E/u))'^  is  a  constant  depending  upon  the  gas.  In  terms  of  the  laser  power 

2  2  I  x\  2  ^  f  x\  2 

required  for  threshold,  P  «  R  E  ,  so  that  P  «  N  R  u  . 

The  exact  dependence  of  the  non-linear  or  multiphoton  initiation 
threshold  on  experimental  parameters  varies  with  the  theoretical  description 
of  the  ionization  process  which  is  most  appropriate  to  the  gas  under  consideration 
Thus,  the  ionization  may  occur  due  to  tunneling  of  an  electron  out  of  the  atom 
through  the  potential  barrier  created  by  the  intense  applied  field. Neglecting 
the  fact  that  the  electric  field  is  alternating,  the  frequency  of  field  Induced 
ionization  events  per  atom  is  given  approximately^  by 
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7  3/2  -1 

a  ■  exp  (-6.8  X  10  1  E  ) 

where  v  is  the  effective  frequency  at  which  electrons  approach  the  potential 
®  16  - 1 

barrier  and  is  roughly  10  sec  ,  I  is  the  ionization  potential  of  the  atom 

and  E  is  the  electric  field  strength.  For  typical  laser  breakdown  conditions 

7  3 

I  ■  15.8  eV  and  E  ■  10  V/cm  so  that  a  ■  10  cm  /sec.  Therefore  field 

g 

strengths  of  the  order  of  3  x  10  V/cm  ifould  be  required  to  produce  significant 
ionization.  This  crude  calculation  would  seem  to  indicate  that  a  multiphoton 
process  is  involved. 

The  ionization  coefficient  per  atom  according  to  the  multiphoton  theory 

is  expected  to  be  proportional  to  (F)’”  where  F  is  the  laser  photon  flux  and 

m  (1/nv).  For  argon  and  a  ruby  laser  m  •  9  >  so  that  the  requirement  that 

3*3 

2:Mt^t  1  means  that  the  threshold  condition  would  require  that  E  MiC^  be  a 
constant  or  that  the  threshold  power  obey  the  relation  Pa  (MAt)  *  R  .  If 
the  gas  ionized  by  the  multiphoton  process  were  mercury  (I  «  10.4  eV),  we 
would  expect  Pa  (MAt)  and  that  the  power  required  would  be  very  much 

lower.  According  to  these  relations  the  power  required  to  meet  this  threshold 
Is  idition  should  be  insensitive  to  changes  in  the  gas  density  and  should  increase 
with  'he  square  of  the  focal  spot  size.  This  means  that  the  multiphoton  process 
threshold  condition  will  become  more  difficult  to  meet  relative  to  the  avalanche 
threshold  condition  as  the  focal  radius  is  increased.  This  should  lead  to  a 
greater  spread  in  the  breakdown  time. 

The  difference  in  behavior  of  the  two  ionization  processes  is  most 
pronounced  when  the  integrated  laser  energy  is  fixed  and  the  spatial  distribution 
of  the  focus  is  varied,  either  intentionally  or  by  accident  as  when  the  laser 
mode  pattern  changes.  Thus  a  20^  increase  in  the  peak  laser  power  at  any  point 
in  space  will  increase  the  ionization  rate  due  to  the  multiphoton  process  by 
factors  of  from  3  for  mercury  to  ^,1  for  argon  whereas  the  ionization  frequency 
due  to  the  avalanche  process  increases  only  by  a  factor  of  about  1.6.  This  high 
sensitivity  of  the  multiphoton,  or  other  non-linear,  ionization  to  the  intensity 
of  the  laser  power  may  account  for  some  of  the  observed  fluctuation  in  the 
breakdown  time.  If  this  effect  is  important  at  a  fixed  laser  pulse  energy,  the 
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shorter  breakdown  times  should  correspond  to  those  pulses  in  which  the  laser 

power  distribution  is  more  sharply  peaked  in  the  focal  region  whereas  the 

longest  breakdown  times  would  correspond  to  pulses  in  which  the  spatial  variation 

in  laser  power  was  relatively  uniform.  Experimental  data  may  be  available  for 

3 

checking  this  hypothesis.  An  alternative  explanation  of  the  observed  spread 
in  breakdown  times  is  given  in  the  next  section. 


3)  Statistics  of  Breakdown 

In  this  section  we  consider  more  explicitly  the  effects  of  a  non-linear 
initiation  mechanism  on  the  statistics  of  breakdown.  Our  main  point  is  that 
the  theory  is  essentially  no  different  from  that  applicable  in  the  case  of 
Townsend  breakdown  when  the  initiating  electron  is  produced  by  a  constant 
external  source  but  with  the  condition  that  the  probability  of  producing  an 
avalanche  is  a  function  of  the  field  strength  or  in  the  case  in  which  the 
initiating  electrons  are  produced  by  field  emission  from  a  cathode.  These 
problems  have  been  considered  previously  by  Phelpo  and  Berg^  and  others.  This 
theory  is  especially  simple  when  the  time  required  for  the  avalanche  to  grow 
to  the  value  required  for  observation  of  breakdown  is  negligible.  In  this  case 
the  probability  of  breakdown  G(t)  after  a  time  t  is  given  by 

in  G(t)  -  -  I  K(C)  dt  (4) 

where  k(t)  is  the  frequency  of  production  of  electrons  capable  of  initiating 
the  avalanche.  According  to  the  laser  breakdown  model  being  discussed  we 
expect  \.(t)  to  be  proportional  to  ME^'"(t),  l.e.,  a  very  rapidly  varying  function 
of  time.  The  probability  that  breakdown  will  be  observed  at  a  time  between 
t  and  t  +  dt  is  found  by  differentiating  Eq.  (4)  to  get 

dG  -  -  \(t)  G(t)  dt  ( 5) 
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This  is  actually  the  basic  equation  for  C(t).  As  shovm  by  Phelps  and  Berg,^ 

Eq.  (5)  can  be  used  to  evaluate  \(t).  Thus, 

-  -  gTT)  f  <*') 

The  behavior  of  the  various  functions  defined  above  is  shown  in 
Fig.  1.  Here  we  have  shown,  in  Fig.  la,  the  time  dependence  of  the  laser 
power,  P(t),  from  an  oscilloscope  photograph  taken  by  Waynant.^  Also  shown  in 
Fig.  la  are  the  relative  values  of  \(t)  as  calculated  for  multiphoton 
ionization  of  helium  (m  ■  11)  and  of  mercury  (m  ■  6).  Figure  lb  shows  curves 
of  G(t)  as  calculated  from  the  X.(t)  for  m  ■  11  with  the  magnitude  of  \(t) 
adjusted  so  as  to  produce  an  average  number  of  electrons  per  pulse,  N,  of  1, 

5,  18,  or  2050.  Here  N  G(o»)  \.(t)dt.  The  corresponding  values  of 

dG/dt  are  shown  in  Fig.  Ic.  The  most  striking  feature  of  these  curves  is  the 
relatively  small  region  of  time,  e.g.,  5*10  nsec  halfwidth,  over  which  the 
initiating  electrons  are  produced  with  an  appreciable  probability.  This  is 
consistent  with  an  apparent  upper  limit  to  the  observed  breakdown  times. 

A  second  feature  of  these  curves  is  the  small  but  definite  shift  of  the  curve 
of  the  probability  of  initiation.  Fig.  Ic,  to  shorter  times  as  the  peak  laser 
power  is  increased.  With  m  ■  11  the  increase  in  N  from  1  to  5.18  and  2050 
represents  an  increase  in  peak  laser  power  by  factors  of  1.16  and  2.0, 
respectively.  If  the  time  required  for  avalanche  growth  were  negligible  then 
the  factor  of  two  increase  in  laser  power  would  result  in  breakdown  at  times 
of  the  order  of  one  half  the  time  required  to  reach  peak  power. 

If  our  assumption  that  the  time  required  for  avalanche  growth  is 

negligible  and  if  fluctuations  in  the  spatial  distribution  of  the  laser  intensity 

can  be  neglected,  then  a  plot  of  the  number  of  breakdown  events  per  unit  time 

3 

should  resemble  Fig.  Ic.  The  data  obtained  for  argon  by  Waynant  show  only  at 
the  higher  laser  energies  are  the  distributions  similar  to  Fig.  Ic.  The  data 
obtained  at  energies  nearer  threshold  yields  distributions  which  are 
considerably  more  spread  out  in  time  and  for  which  the  peak  occurs  much  later 
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in  time.  The  shift  to  later  times  at  lower  laser  energies  is  characteristic 
of  longer  avalanche  growth  constants  and  suggests  strongly  that  the  avalanche 
growth  time  is  not  negligible.  If  the  laser  power  were  constant  during  the 
growth,  then  the  finite  growth  time  could  be  taken  into  account  at  a  fixed  laser 
energy  simply  by  shifting  the  time  scale  in  Figs,  lb  and  Ic  by  an  appropriate 
amount.  With  a  time  varying  laser  power  there  will  be  a  tendency  for  the  cases 
in  which  the  Initiating  electron  occurs  later  to  have  a  lower  average  growth 
constant.  However,  at  present  this  effect  seems  too  small  to  produce  the  observed 
spread.  Obviously,  further  investigation  of  this  area  is  desirable. 


d)  Conclusions 

The  arguments  presented  in  this  report  indicate  that  the  non-linear 
mechanism  which  probably  provides  the  initiating  electron  for  laser  breakdown 
is  multiphoton  absorption.  This  process  produces  electrons  only  over  a 
relatively  narrow  time  interval  and  so  its  effect  is  consistent  with  the 
tendency  of  the  breakdovm  time  to  show  an  upper  limit  when  plotted  as  a  function 
of  laser  power.  The  time  required  for  avalanche  growth  appears  to  become 
appreciable  as  the  laser  energy  is  reduced  toward  threshold  in  argon  at 
1  atmosphere. 
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APPENDIX  B 


A  RANDOM  WALK  PROBLEM  IN  LASER  PHYSICS 
D.  P.  Gaver,  J.  A.  Marshall 


Introduction 

This  report  describes  some  preliminary  work  on  probability  models 
in  physics,  carried  out  at  the  request  of  Dr.  A.  Phelps.  The  results  and 
methods  discussed  provide  a  convenient  starting  point  for  checking  certain 
calculations  and  models  used  by  physicists;  in  particular,  some  easily 
manipulated  "continuous"  approximations  may  be  appraised  for  numerical 
accuracy. 

The  write-up  to  follow  is  in  the  nature  of  a  progress  report. 

More  work  along  lines  similar  to  that  described  here,  can  be  conducted. 

A  First  Model 

A  particle  (electron)  executes  a  random  walk  in  continuous  time 
over  a  denumerable  set  of  (energy)  states  i,  2,  3,  ...  .  Steps  in  the 
walk  (energy  changes)  occur  at  random  times  because  of  interaction  with 
particles  of  another  class  (photons);  the  steps  are  either  up  by  one  unit, 
l.e.  from  state  i  to  i+1,  or  down  by  one  unit,  l.e.  from  i  to  i-1,  except 
that  when  the  lowest  state,  1,  is  reached  no  further  down-jumps  are  possible. 
Furthermore,  when  some  sufficiently  high  state,  H,  is  reached  the  particle 
Jumps  instantaneously  back  to  state  1,  and  the  walk  begins  again.  In  fact, 
when  the  Jump  from  H  to  1  occurs,  another  particle  is  born  and  occupies 
state  1  at  birth,  thereafter  walking  independently  of  the  first.  This 
process  then  continues  indef  ini  ’ y . 
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In  this  report  we  study  the  probability  distribution  of  the 


random  variable  described  as  the  time  of  a  particle's  first  passage 

to  state  H,  having  started  at  state  1.  In  particular  we  compare  the  mean, 

or  expected,  first  passage  times,  E(T,„),  as  the  latter  are  computed  using 
"  In 

various  hypotheses  ab^ut  jump  probabilities.  We  also  discuss  the  problem 
of  continuous  approximation  of  the  path  of  the  particle  through  energy 
space . 


In  order  to  obtain  explicit  expressions  for  E(T,„)  we  make  the 
specific  assumption  that  the  particle  walk  is  a  time -homogeneous  birth- 
and-death  process,  i.e.  that  the  probability  that  a  particle,  currently 
occupying  (energy)  state  i  (i  >  1)  Jumps  to  state  i+l  in  the  time  dt  is 
X^dt  ♦  o(dt);  the  corresponding  down- jump  probability  is  u^dt  +  o(dt), 
except  that  =  0;  importantly,  too,  all  jump  probabilities  are  influenced 
only  by  present  state,  and  not  by  any  previous  history.  This  makes  the 
walk  a  simple  Markov  process.  Both  and  are  positive  and  independent 
of  the  time  during  which  the  process  has  gone  on.  An  Introduction  to  such 
processes  will  be  found  in  Feller  [2].  The  above  assumptions  imply  that 
the  distribution  of  the  sojourn  time,  S^,  of  a  particle  in  a  particular 
state,  i,  i.e.  the  uninterrupted  time  from  an  entrance  into  i  until  the 
next  exit  therefrom,  is  distributed  exponentially 


-(X.+p  )t 

Prob[S^  <t]  =  l-  e  ^  ^ 


(1) 


If,  for  some  reason,  the  exponential  distribution  is  physically  implausibl'* 
the  distribution  (1)  can  actually  be  made  of  arbitrary  form,  leading  us  to 
the  class  of  semi-Markov  rather  than  Markov  processes.  The  methods  to  be 
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described  then  need  only  be  altered  trivially  to  obtain  explicit  solutions, 
provided  there  is  independence  between  the  successive  sojourns.  Independence 
would  not  occur,  for  example,  if  the  event  that  a  particle  had  Just  entered 
state  i  from  below,  i.e.  from  i-1,  were  to  shorten  the  time  of  sojourn  of 
the  particle  in  S^,  or  perhaps  alter  the  probability  with  which  the  particle 
Jumps  from  1  to  1+1  or  to  i-1.  Such  effects  might  derive  from  considerations 
of  velocity  as  well  as  position;  we  rule  them  out. 

Under  our  assumptions,  a  convenient  way  of  computing  the  Laplace- 

"®Tih 

Stleltjes  transform  E(e  ),  the  expected  first  passage  time,  E(T,„),  and 
higher  momenta  as  well  has  been  given  in  Gaver  [3];  formulas  for  have 

been  given  previously  by  others.  We  will  not  utilize  the  L.-S.  transform 
here,  although  by  employing  this  function  it  should  be  possible  to  prove 
■certain  probability  limit  theorems  of  use  when  the  second  stage  of  the 
problem  --  involving  particle  branching  when  Jumps  from  high  energy  states 
to  state  1  —  la  attacked.  From  [  3  ]:  let  denote  the  first  passage 
time  from  state  i  to  i+1  (i.e.,  the  elapsed  time  from  first  entrance 
into  state  i  to  first  entrance  into  state  1+1).  Then  .{' 


a^(8)  =  E(e 


f 


since  initially 

aj^(e) 


X 


i 


8+X^+U 


i-1 


(2) 
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(3) 


we  can  succeaalvely  substitute  to  obtain  a^(s)  for  any  i.  Finally,  because 


Tih  -  +  Ag  +  . . .  4  A^_j^ 

with  A  a  sequence  of  independent  random  variables  we  find 


(4) 


-sT 


E(e  ^«)  = 


H-1 


(5) 


For  the  expectation  or  mean  of  T,  we  use  the  fact  that 


E(Tih)  =  ;  E(A_^) 


H-1 

V 

L. 

i=l 


and  the  recursion  formula,  derivable  from  (3)  by  evaluating  its  derivative 
with  respect  to  a  at  s  =  0, 


E(Ai)  =  ^  E(A^_j^)), 


(6) 


which  in  turn  gives 


E(A^) 


^  ^  ^  1  ^  ^i^i-1  1 _  ^ 

\  \  ^i-1  ^i-2 


(7) 


It  is  expression  (7)  that  will  be  particularly  discussed  in  what  follows. 
Note  that  the  expression  (6)  can  be  easily  evaluated  by  computer,  for  it 
is  a  simple  recursion  formula.  This  has  been  carried  out  in  certain  cases, 
believed  to  be  of  interest  in  applications;  some  examples  will  follow. 

First  we  present  some  special  cases  that  may  be  of  use. 


Special  Case  I:  X.^=X,  (i=l,2,...;  but  |ij^=0). 

In  this  walk  the  probabilities  of  Jumps  and  the  Jamp  rates  (hen. e 
sojourn  times)  do  not  depend  upon  the  states.  From  (6)  or  (7)  ve  easily  get 


E(A^) 


1*  +  ^ 


and  therefrom 


H  =  X,. 


_  H(H-l) 
2X 


(x-n)‘ 


[1  -  (^)  ]  for  X  ^  n 


For  large  H  we  can  put,  for  the  case  X  4  u, 

E(Ti„)  -  —^2 

'■'*  *•  “  (X..M) 

H-1 

*■  (>-ur 

where  the  approximation  Improves  as  H  Increases. 


for  X  =  u. 


for  u  =  X 


for  ^  >  X 


Special  Case  II:  Xj  =  IX,  =  lu  (1=1,2, ...j  but  =0), 

In  this  model  activity  Increases  with  energy.  Again  using  (6) 
or  (7)  we  e7 1 

■  I  i 


If  X  =  u  then  the  above  series  becomes  the  harmonic,  and 


ECAj) 


ilal. 


Therefore 
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(12) 


e(Tih) 


H-1 

r 

L 

1=1 


E(A 


i>~X  / 


H-1 

In  X  dx 


=  i((H-l)  In  (H-1)  -  {H-2)] 

where  the  approximation  Improves  as  H  Increases.  We  will  not  stop  to  write 
down  the  approximations  obtainable  for  X  ^  n;  these  could  be  gotten  by 
replacing  the  sums  by  Integrals,  perhaps  in  a  refined  form  by  use  of  the 
Euler-Maclaurln  sum  formula. 

The  above  models  have  taken  no  account  of  a  necessary  physical 
restriction  on  the  sequences  imposed  by  "detailed  balance". 

This  principle  requires  that  the  transition  rates  for  cases  I  and  II  be 
modified: 


Case  I':  X.  =  X,  |i.  =  ^  —r— 


(13) 


Case  II':  X^  =  IX, 


(i-l)X 


(14) 


Notice  that  for  lerge  i  we  are  very  nearly  in  the  cases  I  and  II,  with 
X  =  u.  However  inclusion  of  detailed  balance  has  the  effect  of  slightly 
decreasing  compared  to  its  corresponding  value  in  cases  I  and  II,  and 
hence  decreasing  E(Tj^j^). 

Computations  were  carried  out  for  cases  I’  and  II'  in  order  to 
compare  the  results  of  the  continuous  model  with  those  jf  the  discrete 
model  developed  above.  Using  (13)  and  (7)  we  find 


E(A^) 


H  tvi 


^^2  . 


\/n  +  sji] 


(15) 
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for  case  I'  while  (l4)  and  (7)  give  l/i  times  this  result  for  case  II'. 

Values  of  the  quantity  F(T  )  were  calculated  for  all  values  of  H  from 

In 

2  to  31  for  each  of  the  two  cases  (with  X  =  1). 

The  corresponding  results  from  the  continuous  model  are 
Case  I':  e(T^jj)  =  (H-l)^/3,  and 

Case  II'.  E(Tj^jj)  =  2(H-1)/3  • 


Table  1  gives  the  results  of  the  computations  for  the  discrete 
model  and  compares  these  with  continuous  model,  for  the  two  cases.  The 
agreement  between  the  two  models  is  seen  to  improve  as  H  increases.  Notice 
that  the  discrete  model  always  gives  the  longer  first  passage  times.  As 
a  check  on  this,  consider  approximating  the  discrete  model  by  another  one 
with  smaller  step  size. 

If  we  want  to  approximate  the  behavior  of  the  particle  ever 

the  original  state  space  by  a  particle  "walking"  on  a  new  state  space  which 

has  smaller  step  size,  we  must  find  the  appropriate  values  for  the  new 

parameters.  If  we  take  the  original  step  size  as  1  and  the  original 

parameters  as  and  and  the  new  step  size  as  A  with  parameters 

X' ,  and  uJi;  then  quantum  mechanics  suggests  for 

1  2 

Case  I':  X|,  =  (i)  X,  emd 


S' 


p 


while  for 


Case  II':  X^ , 


1' 


i'(I)X,  and 
XJ,(1'-1)  vA_^. 
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Table  1 

Values  of  E(T.  „)  for  the  Discrete  and  Continuous  Models 
and  the  Ratio  of  These  Quantities  for  Cases  I'  and  II' 

with  X  =  1. 

Case  I'  Case  II' 


H 

Discrete 

Continuous 

d/c 

Discrete 

Continuous 

D/C 

2 

1.000 

.33 

3.00 

1.000 

.67 

1.50 

3 

2.707 

1.33 

2.03 

1.854 

1.33 

1.39 

4 

5.101 

3.00 

1.70 

2.652 

2.00 

1.33 

5 

8.174 

5.33 

1.53 

3.420 

2.67 

1.28 

6 

11.923 

8.33 

1.43 

4.170 

3.33 

1.25 

7 

16.345 

12.00 

1.36 

4.907 

4.00 

1.23 

8 

21.439 

16.33 

1.31 

5.634 

4.67 

1.21 

9 

27.204 

21.33 

1.28 

6.355 

5.33 

1.19 

10 

33.639 

27.00 

1.25 

7.070 

6.00 

1.18 

11 

40. 744 

33.33 

1.22 

7.780 

6.67 

1.17 

12 

48. 519 

40.33 

1.20 

8.487 

7.33 

1.16 

13 

56.962 

48.00 

1.19 

9.191 

8.00 

1.15 

14 

66.074 

56.33 

1.17 

9.892 

8.67 

1.14 

15 

75.855 

65.33 

1.16 

10. 590 

9.33 

1.13 

16 

86  .  304 

75.00 

1.15 

11.287 

10.00 

1.13 

17 

97. 422 

85.33 

1.14 

11.982 

10.67 

1.12 

18 

109.207 

96.33 

1.13 

12.675 

11.33 

1.12 

19 

121.660 

108.00 

1.13 

13.367 

12.00 

1.11 

20 

134. 781 

120.33 

1.12 

14. 058 

12.67 

1.11 

21 

148. 570 

133.33 

1.11 

14. 747 

13.33 

1.11 

22 

163.027 

147.00 

1.11 

15.435 

14.00 

1.10 

23 

178.151 

161.33 

1.10 

16.123 

14.67 

1.10 

24 

193.9**3 

176.33 

1.10 

16.809 

15.33 

1.10 

25 

210. 402 

192.00 

1.10 

17.495 

16.00 

1.09 

26 

227. 529 

208. 33 

1.09 

18.180 

16.67 

1.09 

27 

245.323 

225.33 

1.09 

18.865 

17.33 

1.09 

28 

263.785 

243.00 

1.09 

19.548 

18.00 

1.09 

29 

282.914 

261.33 

1.08 

20.232 

18.67 

1.08 

30 

302.716 

280.33 

1.08 

20.91** 

19.33 

1.08 

31 

323.173 

300.00 

1.08 

21. 596 

20.00 

1.08 
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Now  consider  the  situation  represented  in  Figure  1  below. 


I  2  3  4  . 

- - - - - - - ‘  i 

•  X 

1  ^  30  40 

Fig.  1 

The  top  line  represents  the  state  space  with  step  size  1  and  H  =  4.  If  we 

take  X  =  1  we  find  from  Table  1  that  E(T,,  )  =  5.101  for  Case  I'.  The  bottom 

line  represents  a  state  space  with  A  =  l/lO.  It  is  believed  that  sliding 

the  first  state  to  the  left  as  shown  is  consistent  with  the  approach  taken 

in  the  usual  continuous  model.  For  the  new  state  space  we  have  (since  X=l) 

=  100.  The  first  passage  time  which  we  want  to  compare  is  that  from 

1  to  30  on  the  new  space.  From  Table  1  we  find  E(Tj^^qj  =  302.710/100  = 

3.027.  This  is  less  than  the  value  for  step  size  1  which  is  consistent  with 

the  fact  that  the  continuous  model  gives  shorter  firs^  passage  times, 

and  it  is  quite  close  to  the  value  of  3.000  given  by  the  continuous  model. 

The  ''.bova  calculations  for  Case  II'  give  =  2.652  for  step  size  1 

/ 

and  E(Tj^^q)  =  20.914/10  =  2.091  for  step  size  l/lO.  Again  the  change  is 
in  the  expected  direction  and  the  continuous  model  gives  2.000. 

It  may  be  of  interest  to  derive  an  approximate  closed  form 
expression  for  E(T^jj)  in  Cases  I'  and  II'  with  X  =  1  and  A  =  1.  In  other 
words  we  would  like  a  closed  expression  for 
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By  the  Buler-Maclaurin  sum 


H-1 

Y 

1=1 

where  E(A^)  la  given  by  (15)  with  X,  =  1. 
formula  (Cramer  [1],  p.  123) 


1  1 
V  //2 
z. 

V=1 


,-V2 


=  Y dx  >  I  +  I  1^/^  -  I  Y Pi(x)x""'^  dx 


where  P^(x)  =  (xj  -  x  +  ^  and  [x]  Is  the  largest  Integer  not  larger  than  x. 
Now  P^(x)  <  ^  ao 


(x)x'^^^dx  S  A 


and  therefore 


1  1 

V  1/2 

v=l 


1 

r  1/2  ^  11  ,1/2  1/,V2  T  X 

>  /  X  '  dx  +  2  +  2  ^  -  2^^ 


-  S  l2/3  .  2  ^  1  .  2  i2/3 
3  3  3 


Thus  for  Case  I'  we  have 

H-1 


E(T,,)  .  Y  e(a,)  >Y  I 


1=1 


1=1 


which  Is  to  be  compared  with  the  result  (H-1)  /3  obtained  from  the 
continuous  model.  For  Case  II', 


1=1 


1 =1 
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which  is  exactly  the  result  obtained  from  the  continuous  model. 

“nils  approach  could  perhaps  be  used  with  a  better  estimate  of 
P-(x)  to  give  a  more  accurate  expression  for  E(T, „)  than  is  obtained  with 
the  continuous  model. 

Continuous  Approximations 

It  is  well  luiown  that  the  path  of  a  classical  random  walk, 
itself  a  step  function,  may  often  be  usefully  approximated  by  a  continuous 
function.  Under  many  circumstances,  this  results  in  simplicity,  and  such 
an  approach  hr  3  been  used  to  discuss  the  present  problem.  We  shall 
explain  the  approximation  method  commonly  used  in  terras  of  our  Case  I 
above.  In  terms  of  that  example  we  can  proceed  heuristically  in  a  natural 
way. 

Rather  than  pass  immediately  to  the  continuous  function  approxima¬ 
tion,  let  us  suppose  that  our  original  process  steps  are  of  unit  size 
(as  they  are),  and  that  we  wish  to  pass  to  an  equivalent  process  acting 
over  energy  states  of  step  size  ^  .  The  idea  of  such  equivalence  is  as 
follows;  suppose  our  original  particle  begins  its  random  walk  on  the 
original  energy  states,  l.e.  on  those  differing  by  unity;  here  we  disregard 
the  boundaries.  Then  at  time  t  after  the  walk  begins  its  position  will 
have  a  certain  probability  distribution:  the  expected  distance  from  the 
starting  point  and  the  variance  of  the  resulting  distribution  will  be 
E(Dj^)  =  (X  -  u)t 

and  ( l6 ) 

Var(D^)  =  (X  ^  M)t 
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and  if  X.t  and  nt  are  at  all  large  then  the  distribution  of  tends  to 
the  Gaussian  described  by  the  above  parameters.  If,  on  the  A  scale,  we 
start  a  walk  from  the  same  place,  utilizing  rate  parameters  X  and  u 
respectively,  we  find  for  ,  the  displacement  on  the  A  scale 


and 


E(D^)  =  (X  -  u)At 

Var(D  )  =  (X  +  m)  • 

A 


(17) 


(18) 


Again  tends  to  the  Gaussian  form  with  the  same  parameters.  We  say  that 

the  -process  is  equivalent  to  the  Dj^-process  if  the  above  two  limiting 

(kiussian  distributions  are  the  same,  or  if 
E(D,)  =  E(D^  ) 

^  A 

Var(Dj^)  =  Var(D^  ), 

which  implies  that 

X  -  u  =  (X  -  u)A 

X  +  u  =  (X  +  u)A  , 


(19) 


so  we  can  solve  for  X  and  u  to  find 


X 


it 


A  A2 


U 


2^  A 


■ 


(20) 


It  is  intuitively  appealing  that  any  probability  question  asked  about 
equivalent  processes  should  provide  approximately  equal  answers  -- 
particularly  if  both  processes  must  operate  "for  a  long  time",  for  this 
allows  the  tendency  towards  the  Gaussian  limiting  form  to  occur.  Thus  we 
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may  expect  that  we  can  utilize  processes  on  finer  and  finer  step  sizes  to 
calculate  first-passage  times  to  a  "high  state"  H,  with  better  approximation 
as  H  increases.  We  are  interested  in  this  possibility  because  it  is 
conventional,  and  often  convenient,  to  pass  to  infinitesimal  st'p  sizes 
(let  A — ^  0)>  to  replace  the  random  walk  by  "diffusion"  or 

Brownian  motion.  We  are  interested  in  the  error  in  such  an  approximate 
procedure.  Refer  to  the  following  figure,  illustrating  the  original  walk, 
and  the  new  one: 


W 


C  I  Z 


Fig.  2 

Let  us  examine  our  Special  Case  I,  translating  to  the  new  scale 
with  step  size  ^ .  Take  ^  to  be  an  integer,  so  that  we  have  an  integral 
number  of  ^-steps  replacing  one  l-step.  Further  let  both  walks  start  at 
the  same  place  on  the  original  scale:  thus  the  new  walk  will  start  at 
i  =  A  for  i  on  the  new  scale,  and  will  end  at  hA  Then  from  (6) 
and  (20)  we  get,  for  1  on  the  new  scale. 
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1 


^  1  -(^)' 

X 


E(A^) 


^  1  -  ^ 


1 

X 


M  4  >■ 


(21) 


M  =*  X 


and 


e(t^j^IA)  = 


hA‘^-1  _  .  _  A 

I  fi  -  &]-  ■  & 

X-M  X  1  X-n  ^  ^  X 


-11- 


-  (h-i)A'^ 

(-) 


(h-i)(h.i-A) 

2l 


X  ^  M 

X  =  ^  (22) 


Now  the  appropriate  scaling  is  given  by  (19)  and  (20).  When  we  substitute 
for  X  and  u  in  terms  of  the  original  X  and  n  a  complicated  expression 
results  for  general  A*  However,  when  A — ^  0  a  relatively  simple  limit 
results; 


E(T^jj|0) 


H-1 

X-p 


+ 


1 

2 


X-fp 

(x-p)^ 


-2(^)H  -2(^) 

g  ^Xfp'  ^X+p' 


H^-1 

2X 


^  =f 


(23) 

X  =  p  . 


One  can  see  that  E(Tj^jj|  A)  increases  for  decreasing  A  when  X  =  p;  this 
shows  that  for  this  case  the  continuous  approximation  over-estimates 
time  of  first  passage  to  H.  Although  the  comparison  for  different  X  and 
p  caji  be  made  using  (22),  we  have  not  done  so  here. 
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